Interplay between Superconductivity and Antiferromagnetism in a Multi-layered 

System 
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Based on a microscopic model, we study the interplay between superconductivity and antifer- 
romagnetism in a multi-layered system, where two superconductors are separated by an antiferro- 
magnetic region. Within a self-consistent mean-field theory, this system is solved numerically. We 
find that the antiferromagnetism in the middle layers profoundly affects the supercurrent flowing 
across the junction, while the phase difference across the junction influences the development of an- 
tiferromagnetism in the middle layers. This study may not only shed new light on material design 
and material engineering, but also bring important insights to building Josephson-junction-based 
quantum devices, such as SQUID and superconducting qubit. 

PACS numbers: 74.50.+r, 75.70.-i, 74.81.-g 



I. INTRODUCTION 

As a remarkable aspect of high-T c superconductivity, 
its unique characteristics may result from the competi- 
tion between more than one type of order parameter r^— 
Historically, the typical antagonistic relationship between 
superconductivity and magnetism has led researchers to 
avoid using magnetic elements, such as iron, as potential 
building blocks of superconducting (SC) material. How- 
ever, it is now well accepted that the unconventional 
superconductivity emerging upon doping is closely re- 
lated to the antiferromagnetism (AFM) in parent com- 
pounds 4-7 . Recent research theme in high- J), cupratc 
community centers on how to establish the connection 
between antifcrromagnctic (AF) and d-wave SC (DSC) 
orderings^, i.e., whether they compete with each other 
or coexist microscopically. This issue is a subject of cur- 
rent discussions. Depending on material details, some 
compounds show the coexistence of DSC and AFM2 
while others exhibit microscopic separation of these two 
phases^. Engineered heterogeneous systems and multi- 
layered high-T c cuprates offer a unique setting to study 
the interplay between DSC and AFM. In these systems, 
the disorder effect can be minimized significantly with 
atomically smooth interfaces. Experimentally, no mix- 
ing of DSC and AFM was reported in the heterostruc- 
ture artificially grown by stacking integer number of 
Lai.85Sr .i5CuC>4 and La2Cu04 layers 1 ^. Meanwhile mi- 
croscopic evidence for the uniform mixed phase of AFM 
and DSC in outer Cu02 planes was reported^ on a Hg- 
based five-layered cuprate. Theoretically, Demler et alr^ 
have studied the proximity effect and Josephson coupling 
in the SO (5) theory of high- T c superconductors. Depend- 
ing on the thickness, the middle antiferromagnetic region 
could behave like a superconductor, metal or insulator. 
On the other hand, it was shown in the perturbation the- 
ory 1 ^ that the spin exchange coupling in the insulating 
AF layer can allow the tunneling of Cooper pairs. 

In this article, we study microscopically a multi-layered 
system with two superconductors separated by AF lay- 
ers. Through self-consistent mean-field theory and nu- 



merical diagonalization, we are able to solve the Hamil- 
tonian, which allows us to study the interplay between 
DSC and AFM in detail. Varying some external param- 
eters, such as the SC phase difference across the system 
Sip and Coulomb interaction U in the AF layers leads to 
interesting results about the interplay. 




FIG. 1: Schematic drawing of a multi-layered system. The 
middle region consists of AF layers (in X-Y plane) and the 
two sides are SC layers with d-wave pairing symmetry. The 
phases of the SC order parameter (or the pairing potential) 
in the first (far left) and the last (far right) layer are fixed at 
ifL and ipn, respectively. 



II. MODEL AND SETUP 

The model system under consideration is schematically 
presented in Fig. [T] for which the Hamiltonian can be 
written as: 



ijTllTlj fj> y Til , 



in 



where denotes the hopping integral between the near- 
est neighbor sites. For in-plane (X-Y plane) hopping, 
ti j = t, and for inter-layer (Z direction) hopping, tij = 

t±- Ci,a (cJ CT ) is the annihilation (creation) operator of 
electrons on the ith lattice site with spin a (a =j, 4-). 



2 



The quantity rij l(T = c\ a Ci l(T is the number operator on 
the ith. lattice site. Both C/j and are positive. Ui indi- 
cates the on-site repulsive Coulomb interaction, which is 
nonzero only in the middle layers, and Vy describes the 
in-plane nearest neighbor attractive interaction, which is 
nonzero only in the SC layers. We are going to fix the 
phase of the SC order parameter in the first (far left) 
and the last (far right) layer at tpi, and tp^, respectively, 
which can be introduced by a gauge flux. 

When there is no inter-layer coupling, t± = 0, the 
system is decomposed into individual two-dimensional 
(2D) systems. It is known that in a 2D SC layer, the 
in-plane nearest neighbor attractive interaction Vij leads 
to a nonzero energy gap A with d-wave symmetry, and in 
a 2D AF layer the on-site repulsive interaction Ui leads 
to a nonzero (n, 7r) spin density wave (SDW) order. In 



a genuine cupratc compound, there exists the inter-layer 
coupling, which is usually weaker than the in-plane hop- 
ping. Because of the inter-laycr coupling, there arises 
interesting interplay between SC layers and AF layers 
(see Fig. [J). 

The strategy of solving the Hamiltonian is briefly sum- 
marized as follows: Due to the translational symmetry 
in the 2D X-Y plane, the 3D problem can be decom- 
posed into a 2D (X-Y) plus ID (Z) problem. In the X-Y 
plane, k = {k x ,k y ) is a good quantum number, where 

_I/V +1 ... IN — 1 n — — -N -i/V +1 ••• 
hN y — 1. In order to simplify the Hamiltonian, we adopt 
the mean-field approximation, which leads to the follow- 
ing Bogoliubov-dc Gcnncs equation^: 
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where m and n label the layer number; k and k + Q with 
Q = (it, 7r) denote the wave vectors in the first Brillouin 
zone of the X-Y plane, t mn = t ± (S n , m+1 + <5„, m _i) de- 
scribes the nearest neighbor intcr-plane hopping, ^ = 
— 2t(cosk x + cos k y ) — n is the normal metal disper- 
sion relation, [i is the chemical potential, the variables 
U n is equal to U in the AF layers and zero otherwise 
while V n is equal to V in the SC layers and zero oth- 
erwise. The superconducting pairing potential A^ n = 
V n *i> n (cos k x — cos k y )/ '2, where is the superconduct- 
ing pairing wavefunction. Obviously, the superconduct- 
ing pairing potential, which is the SC order parameter 
has the same phase as that of the superconducting pair- 
ing wavefunction. M n = (n n ^) — describes the 
SDW in the n-th layer. The average electron number 
(n„.-|-), {n n ,i) and the SC pairing wavefunction ^„ can 
be determined sclf-consistcntly through iteration over the 
following relation: 



2 ^ (cosfcz — cos k y ) f w _ 1 f3E, 
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u°i v"'* -u? - v?*~ , and/(£ Q ) = l/(l + eJ 3E «) 
is the Fermi distribution. Through this simplification, 
the multi-layered system is then solved numerically We 
are especially interested in the interplay between the DSC 
and the AFM, which is characterized by the SDW M n 
and the SC pairing wavefunction We will focus 

on zero temperature T = 0. The size of X-Y plane is 
N x = N y = 40. The total layer number is N z = 22. The 
two middle layers are AF layers (U n — U, V n = 0, for 
n = 11, 12), and the rest are SC layers (U n = 0, V n = V, 
for n = 1, 2, • • • , 10, and n = 13, • • • , 22). For simplicity, 
we choose t = 1, t± = 0.1, V = 1.5, and fi = 0. We will 
vary the control parameters, such as U and Sip = ipL — ^PR, 
to study the competition between DSC and AFM. When 
the phase difference Sip is introduced, there is a current 
flowing across the system and its expression is given by: 
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f(E a ) , (3b) w here /± = f(±E a ). 



f(-E a ) , (3c) HI. INFLUENCE OF ANTIFERROMAGNETISM 

ON SUPERCONDUCTIVITY 



where N = N x x N y is the sites number in a 2D plane. 
k samples half of the first Brillouin zone. g(u,v) = 



In the model introduced above, the change of 
the Coulomb interaction U in the middle layers 
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can drive a mctal-Mott insulator phase transition at 
a certain value U c . Correspondingly, the multi- 
layered system changes from a superconductor/normal 
metal/superconductor (SNS) weak linki£ to a su- 
perconductor/insulator/superconductor (SIS) Josephson 
junction^. The studies of SNS and SIS junctions with 
static potential barrier have been well documented^ - — . 
Here we start from a microscopic model and drive the 
middle layers to change from one electronic state into an- 
other by tuning the on-site Coulomb interaction U. We 
are interested in how the SC pairing wavefunction and 
the current across the junction change with the Coulomb 
interaction U. We plot the current as a function of the 
phase difference Sip in Fig. 2. It can be seen that below 
a threshold value of U c « 0.8, the current is a piecewise 
periodic function. In each periodic region, it varies lin- 
early with Sip. This agrees with the result obtained in 
Ref. [IB for a junction with a normal metal. We also find 
that when the Coulomb interaction is larger than U c , the 
current as a function of Sip changes gradually with the 
increase of U. When U is in the range 0.8 < U < 1.5, 
the current shows a shape intermediate between a piece- 
wise linear function and a sinusoidal function. Similar 
results have been observed in Refs^ ' 13 ' 16 when one varies 
the temperature or the thickness of the middle layers in- 
stead of the the Coulomb interaction U . When U « 2.0, 
the current is well approximated by a sinusoidal func- 
tion of Sip (see Fig. 2), which is a typical feature of the 
dc Josephson junction (JJ) current^. The magnitude of 
current decreases rapidly with the increase of U. 

In order to have a better understanding of the influence 
of AFM on DSC, we plot the absolute value and the phase 
of the SC pairing wavefunction in Fig. [3] It can be seen 
that when U is small, there is a nonzero paring poten- 
tial in the middle layers due to the proximity effect, and 
the two superconductors are weakly linked by the nor- 
mal metallic layers. In this weak link regime, the phase 
of the SC pairing wavefunction is almost linear in layer 
number n. This means that when the middle layers are 
in the normal metallic state, it has little influence on the 
SC layers on two sides. However, when U is larger than 
U c ~ 0.8, SDW develops, and the middle layers become 
Mott insulator. Quantum mechanically, Cooper pairs can 
still tunnel through the SDW region when U is not very 
large. We call this intermediate U regime as the tun- 
neling regime. From the weak link regime to tunneling 
regime, the phase of the SC paring wavefunction changes 
gradually from a (nearly) linear function to a step func- 
tion of the layer number (see Fig. [3]). This corresponds 
to the observation that the current shape changes from 
being piecewise linear to sinusoidal. Therefore, when U is 
very large, the current-phase relation becomes identical 
to the famous Josephson rclatio n 17 i 19 . Here for the first 
time, we establish the connection between the profile of 
the SC phase distribution (Fig. [3J and the current-phase 
dependence (Fig. [2J , which we believe is rather intrigu- 
ing. 

When one continues to increase U, the middle layers 



I/10" 3 




0.2 0.4 0.6 0.8 1.0 




FIG. 2: Current across the junction as a function of the phase 
difference Sip of two superconductors for different U. Here 
V = 1.5, t± = 0.1, iV* = N y = 40, N z = 22. Up panel: 
U = 0, 0.1, 0.2, • ■ • , 1.5, and Sip = ~ n. Down panel: 
U = 2.0 and Sip = ~ 2tt. 

becomes more opaque and the superconducting pairing 
wavefunction is completely suppressed in this region. As 
a result, the tunneling of Cooer pairs will be completely 
quenched, and the two comprising superconductors be- 
come decoupled. 

We also would like to mention that when the middle 
layers are in the normal metallic state, the current across 
the junction is a mesoscopic effect^,. The current is pro- 
portional to the phase gradient of the SC pairing wave- 
function. That is, when one increases the number of SC 
layers, the current will decrease and finally vanish. How- 
ever, when the middle layer is in the AF state, the whole 
system becomes a J J. The current of a J J is no longer 
a mesoscopic effect. It is completely determined by the 
phase difference of the SC pairing wavefunction on both 
ends of the system. Therefore, it will not decrease with 
the increase of the number of the SC layers. In this sense 
we say that AFM insulator layers enhance supercurrent. 



IV. INFLUENCE OF SUPERCONDUCTIVITY 
ON ANTIFERROMAGNETISM 

In the above discussion, we find the profound influence 
of the Coulomb interaction U (and hence AFM) of the 
middle layers on the DSC. We now study the back-action 
of DSC on AFM. We will fix the phase difference Sip at 
different values, and see if the SDW in the middle layers 
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FIG. 3: Absolute value (Up panel) and argument (Down 
panel) of superconducting pairing wavefunction 9 n as a func- 
tion of layer number. The phase of SC pairing wavefunction is 
the same as that of the pairing wavefunction. Here V = 1.5, 
t± = 0.1, N x = N v = 40, iV 2 = 22, Sip = 3tt/4. U is fixed at 
different values: (7 = 0, 0.1, 0.2, ■ ■ • , 1.5. 



changes with Sip. We plot in Fig.@]the SDW as a function 
of U for various values of Sip. Clearly, the SDW in the 
middle layers is influenced by Sip in certain range of U. 
When Sip = tt, the critical value U c , after which the SDW 
develops, is much smaller than in the case of Sip = 0, 
though the current is zero in both cases. This is very 
similar to the result in Ref. [TH that the middle AF layers 
will be influenced by the current fed into the junction. 
When the Coulomb interaction is very weak, the SDW 
does not develop for 5<p 7^ n, and the middle layers are in 
the normal metallic state. In contrast, when the Coulomb 
interaction is very strong, e.g. U > 1.2, the AFM is 
robust against the phase difference across the junction 
and strongly suppresses the pairing wavefunction in the 
middle region (see Fig. [3]) . Only when U is small, is the 
emergence of the SDW sensitive to the phase difference. 

When switching off the interlayer coupling (t± = 0), 
the system is decoupled into 2D systems. The DSC on 
two sides will not influence the AFM in the middle layer. 
The middle layer is described by a 2D Hubbard model, 
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FIG. 4: (Color online) SDW in the middle layer (n = 11) 
as a function of the Coulomb interaction for different fixed 



Sip. 



Here U 

= 22, Sip 







1.5, V = 1.5, t± = 0.1, N x = N y 



40, 



(solid), 37r/4 (dashed), and n (dot dashed). 



in which the SDW will develop for an infinitesimal U. In 
the present model (tj_ ^ 0), when Sip = 0, the SDW in 
the middle layers is suppressed due to the inter-layer cou- 
pling, which weakens the perfect nesting at Q = (tt,tt). 
The SDW will not develop until U is larger than a thresh- 
old value U c (see Fig. Hence both t± and 6<p will in- 
fluence the AFM. The results in Fig. @] indicate that the 
AFM shows a crossover behavior for Sip = tt rather than 
a criticality as obtained for Sip = 0. For Sip 6 (0,7r), the 
SDW critical point U c is less than that for Sip = 0. This 
result is in agreement with that in Ref. [l3l 

In order to have a better understanding of the influence 
of the DSC on the AFM in the middle layers, we plot in 
Fig.[5]thc local density of state (LDOS), as given by 



n,t fc+<2,n,t 



k . a 



3/M 



k,a> 



duj' 



(5a) 
05b) 



where oj = E — E a and u>' = E + E a . From Fig. [3 we can 
see a sharp intensity decay of LDOS at Ef when Sip = 0. 
When Sip = tt, the LDOS exhibits a peak structure at Ef- 
Since the development of the SDW is also determined by 
the normal-state DOS at low energies, it explains the 
different SDW critical behavior between the cases with 
Sip = tt and Sip = 0. 



V. CONCLUSION AND DISCUSSION 

We have studied the interplay between the DSC and 
AFM in a multi-layered system. It is revealed that the 
AF layers have profound influence on the DSC. A metal- 
Mott insulator phase transition makes the multi-layer 
system change from a SNS to a SIS junction. We have 
also established the connection between the profile of the 
SC phase distribution and the current-phase dependence. 
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FIG. 5: (Color online) LDOS of the middle layer (n = 11) 
pu,i(E). Here U = 0, V = 1.5, t x = 0.1, N x = N y = 40, 
N z = 22, <Sv? = 0, and tt. 

We find that when the barrier of the middle layers is not 
high enough (metal or semiconductor regime), the mid- 
dle layers cannot maintain a step-function-like phase of 
two superconductors on two sides. Accordingly, the cur- 
rent across the junction will deviate from a sinusoidal 
function of the phase difference. When the barrier of 
the middle layer is high (deep insulator regime), as pre- 
dicted by Josephsoni^ the phase across the junction is 
a step function, and the current is a sinusoidal function 



of the phase difference. Meanwhile, the phase difference 
across the system will dramatically influence the devel- 
opment of SDW in the middle layers. The results from 
our simulations may shed new light on material engineer- 
ing. The DSC and AFM can coexist in a single layer, but 
they also compete with each other. Another important 
insight from the present study is the microscopic model 
of J J. J J has been widely used in a lot of quantum de- 
vices, such as in SQUID for ultra-sensitive magnetic field 
measurement and superconducting qubi1>2i for quantum 
computing. Our study provides a microscopic theory for 
J J, hence has possible applications in building adaptive 
JJ based devices by engineering the middle layers of the 
J J. Finally, with the development of experimental tech- 
niques, it is now possible to synthesize heterogeneous 
systems with layer by layer at atomic scaled. We ex- 
pect that our results can be experimentally verified, as 
some experiments on engineered material have been re- 
portecUi^^. 
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